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We study the role that tachyon fields may play in cosmology as compared to the well-established 
use of minimally coupled scalar fields. We first elaborate on a kind of correspondence existing 
between tachyons and minimally coupled scalar fields; corresponding theories give rise to the same 
cosmological evolution for a particular choice of the initial conditions but not for any other. This 
leads us to study a specific one-parameter family of tachyonic models based on a perfect fiuid mixed 
with a positive cosmological constant. For positive values of the parameter one needs to modify 
Sen's action and use the cr-process of resolution of singularities. The physics described by this model 
is dramatically different and much richer than that of the corresponding scalar field. For particular 
choices of the initial conditions the universe, that does mimick for a long time a de Sitter-like 
expansion, ends up in a finite time in a special type of singularity that we call a big brake. This 
singularity is characterized by an infinite deceleration. 

PACS numbers: 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 

The recently discovered cosmic acceleration Q, |^ I3 has put forward the problem of unravelling the nature of the 
so called dark energy responsible for such phenomenon (for a review see 0, |[5,j ) . The crucial feature of the dark energy 
which ensures an accelerated expansion of the universe is that it breaks the strong energy condition. 

The tachyon field arising in the context of string theory 0, 0, B provides one example of matter which does the 
job. The tachyon has been intensively studied during the last few years also in application to cosmology Ili|-f3^; 
in this case one usually takes Sen's effective action Q for granted and studies its cosmological consequences without 
worrying about the string-theoretical origin of the action itself. We take this attitude in the present paper. 

However, it would be reasonable, before considering concrete tachyon cosmological models, to answer a simple 
question: is the tachyon field of real interest for cosmology? Indeed, it is well-known that for isotropic cosmological 
models, for a given dependence of the cosmological radius on time it is always possible to construct a potential for 
a minimally coupled scalar field (in brief: scalar field) model, which would reproduce this cosmological evolution 
(see e.g. [SJ), provided rather general reasonable conditions are satisfied. Since a similar statement holds also for 
cosmological models with tachyons, one can find some kind of correspondence between minimally coupled scalar field 
models and tachyon ones describing the same cosmological evolution Therefore, a natural question arises: does 
it make sense at all to study tachyon cosmological models in place of the traditional scalar field models? 

In our opinion, the point is that the correspondence between tachyon and scalar field cosmological models is a 
rather limited one and amounts to the existence of "corresponding" solutions of the models, obtained by imposing 
certain special initial conditions. If one moves away from these conditions the dynamics of the tachyon model can 
become more complicated and very different from that of its scalar field cousin. 

In this paper we consider some examples of scalar and tachyon field isotropic cosmological models having coincident 
exact solutions. All the exact solutions considered here actually arise as solutions of some isotropic perfect fluid 
cosmological models. 

Some of the tachyon and scalar field potentials that we consider here are well known in the literature and are 
widely used as quintessential models. We will introduce also a new (at least to our knowledge) tachyon model that 
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is based on the cosmic evolution driven by the mixture of a cosmological constant and of a fluid with equation of 
state p = ke, with — 1 < fc < 1 (e and p denote as usual the energy density and the pressure respectively). The 
corresponding tachyon potential V(T) is represented by a rather cumbersome trigonometric function of the tachyon 
field. When the parameter k is positive the potential becomes imaginary passing through zero for finite values of the 
tachyon field T. Thus, it might seem that these features would kill the model; on the contrary, using the methods of 
the qualitative theory of ordinary differential equations, it is possible to extract from the model a dynamics which is 
perfectly meaningful also for positive fc. We find that this dynamics is quite rich. In particular, when the value of the 
parameter fc is positive, the model possesses two types of trajectories: trajectories describing an eternally expanding 
universe and those hitting a cosmological singularity of a special type that we have chosen to call a big brake. Thus, 
starting with a simple perfect fluid model and trying to reproduce its cosmological evolution in scalar field and 
tachyon models, we arrive to a relatively simple scalar field potential with a correspondingly simple dynamics and to 
a complicated tachyon potential. The latter provides us with a model having a very interesting dynamics giving rise 
to very different cosmological evolutions and opening in turn opportunities for some non-trivial speculations about 
the future of the universe. 

The structure of the paper is as follows. In the second section we study the problem of the correspondence between 
tachyon and scalar field cosmological models and give some examples in the third section. In the fourth section we 
suggest a way to go beyond the limitations of tachyonic models and we make use of this idea in the fifth and sixth 
sections where we study the dynamics and the cosmology of a particular (toy) tachyon model introduced in Sec. 3. 
In the last section we conclude our paper by some speculations on the future evolution of the universe arising from 
the analysis of our model. 



II. THE CORRESPONDENCE 



We consider a flat Friedmann cosmological model ds^ ~ dt^ — a^{t)dl'^ of a universe filled with some perfect fluid 
and suppose that the cosmological evolution a = a(t) is given; the Friedmann equation 



d2 



(1) 



provides the dependence e = e{t) of the energy density of the fluid on the cosmic time (we have put SvrG/S = 1 for 
convenience). Then, the equation for energy conservation 



-3-{e+p) (2) 
a 



fixes the pressure p = p{t); therefore an equation of state p = p{e) can in principle be written describing the unique 
fiuid model compatible with the given cosmic evolution (provided that £ 0). 

Now, let us suppose that the matter content of the universe is modeled by a homogeneous tachyon field T{t) 
described by Sen's Lagrangian density 



L = ^V{T)^/l-t^. (3) 
The energy momentum tensor is diagonal and the corresponding field-dependent energy density and pressure are given 

by 

VI -T2 



p = -T/(r)Vl-T2, (5) 
while the field equation for the tachyon is written as follows 

f 3af V.T 



1 - T2 a V 



= 0. (6) 



One can try and find a potential V{T) so that, for certain suitably chosen initial conditions on the tachyon field, 
the scale factor of the universe is precisely the given a{t). A similar construction can be attempted for a minimally 
coupled scalar field, described by the Lagrangian density 
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A given cosmological evolution a = a{t) can therefore be used to establish a sort of correspondence between the two 
potentials V{T) and U{ip) (see also (1J|). 

Let us see in some detail how this works and what the sought correspondence really means. It is often more practical 
to use the scale factor a to parameterize the cosmic time; this can be done provided d 7^ 0. From Eqs. Q and ||SJ) it 
follows that 



= (8) 



e 

By using Eqs. and Q the latter equation can be rewritten as follows 



where "prime" denotes the derivative with respect to a. Once a{t) and therefore £(a) is specified (see Eq. Q), Eq. 
@ can be integrated to give 



rr^ \ r -e'{x)x 

(an arbitrary additive integration constant is hidden in the unspecified integration limit). By inverting Ea. HlUf) 
a = $^^(r) and making use of the relation 



(£2a6)/ 



the shape of the required tachyon potential V = V{a) = V{^^^{T)) can be found. 

For the minimally coupled scalar field we have similar formulas. The analogs of Eqs. Q and Q for the scalar field 
are 

e=\^'' + U{^), (12) 



l^'-U{^). (13) 



From these we get 



and 



a V e a V 3e 



Integration of Eq. (|14|l gives 

rdx l ~e'{x)x 

and as before U — U{a) — U{F^^{(p)). The formulas above establish a kind of correspondence between the potentials 
U and V in the sense that the cosmologies resulting from such potentials admit the given cosmological evolution 
a = a{t) for suitably chosen initial conditions on the fields. However, it has to be stressed that for arbitrary initial 
conditions on the fields the cosmological evolutions (within the corresponding models) may be drastically different. 
Furthermore, changing the initial conditions, say, for a minimally coupled scalar field theory, one gets different 
cosmological evolutions which, in turn, can be reproduced in entirely different tachyon theories and viceversa. Thus, 
any scalar field potential has a whole (one-parameter) family of corresponding tachyon potentials and the same is true 
the other way around. 
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We can learn something more by expressing the fields T and ip and the potentials V{T) and U{(p) in terms of the 
Hubble variable h — a/ a. Since 

h = -\{e+p) (17) 

one has easily that 
and that 

=s+p = -'^h, U{ip) = h^ + ^h. (19) 

Eqs. p8|l and (|19|l call for the condition e + p > 0, oi equivalently h < 0. For the tachyon model Eq. (|18|l requires 
the additional condition 

h > -^h"^ (20) 

which is equivalent to ask for a negative pressure p < 0. It seems therefore that tachyonic models are more restrictive 
than minimally coupled scalar fields in their possibility to describe cosmological evolutions. We will see in an example 
that this seemingly negative conclusion can be overcome and pleasant surprises may arise. 



III. EXAMPLES 



We now illustrate the above considerations by means of some explicitly solvable examples. In these examples we 
assume, as usual, that the initial moment of time, t = 0, corresponds to the cosmological singularity. All these 
examples are based on simple models of perfect fluids as specified by their equations of state. 

1. The starting point of the first example is a model of universe filled with a perfect fluid with equation of state 
p — ke, with — 1 < fc < 1; the Hubble variable of this model is given by 

A minimally coupled scalar field theory that produces the same evolution h{t) for suitably chosen initial con- 
ditions can be based on one of the following scalar potentials, which arc obtained by the procedure outlined in 
the previous section: 

C;^(^)^2_l_|_^^3^(„). (22) 

the limiting case k ~ 1 gives a minimally coupled massless field. For fc > 1 the potential becomes negative 
reflecting the fact that the velocity of sound is bigger than the speed of light. 

Both choices of sign in the exponent are acceptable and the normalization of the potentials can be chosen 
arbitrarily (the constant tpo)] then, the exact solutions of the field equations providing the prescribed cosmological 
evolution l|21|l are, respectively, the following: 

'P±(t)=T-7^1nt + ^o. (23) 
3v 1 + fc 

The theories are obviously connected by the symmetry operation ip ^ ipo —{ip — ipa), which exchanges the 
potentials and the corresponding solutions. 

To construct a tachyonic field theory we have to restrict our attention to the case p < 0, i.e. fc < 0. Following 
the same procedure we get a tachyon field model, based on the potential 
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There are again two exact solutions of the field equations H21|l . which now do co-exist within the same model: 



T±{t) ^±VT+k t + To. (25) 

This tachyonic model has been studied in |l3l IT^ , where also the correspondence with a minimally coupled 
scalar field theory with exponential potential was noticed. There exists a vast literature on the latter (see e.g. 

MM)- 

2. Now we add a positive cosmological constant to the previous model, i.e. we consider the mixture 

pi = kei, ~1 < k <1, PA = -£A = -A, A > 0. (26) 
This is the same as one fluid with state equation 

p= he - {l + k)A. (27) 
The evolution of the Hubble variable is now given by 



rr 3VA(1 + k)t 
h{t)^VA coth ' . (28) 

The corresponding minimally coupled scalar field theory is based on the following potential: 



There are two exact solutions of the field equations reproducing the given cosmic evolution H28|l : 

(t ) = ± y= In tanh ^-^^ ^ + (^o ■ (30) 

To find the tachyonic model we observe first of all that to have the restriction p < satisfied it is sufficient to 
require that — 1 < fc < 0; this condition is also necessary if one demands that p < during all the stages of the 
cosmological evolution and for any choice of the initial conditions in the tachyonic model. With this condition 
we obtain the following more complicated tachyon potential: 



V{T) = ^ T - (1 + ^) cos^ 



2 

sm 



fVA(l + fc) (T-To) 



-VA(l + fc) (T-To) 



(31) 



Still, the corresponding exact solutions can be found and are given by 



2 . , 3\/A(l + A:)i 

Tit) = ± , arctansinh i '- + Tq. (32 

^ ' 3v/A(l + A:) 2 ^ ^ 

For fc > there are values of the tachyon field for which the potential (I31|l becomes imaginary passing through 
zero. 

When A tends to zero one expects to recover the model studied in the first example. This is indeed the case for 
the tachyonic field. 

As for the scalar potential H29|) the situation is a bit trickier: the correct limit can be obtained only by letting 
the constant ip^ vary with A as follows: 

1 A(l-fc) 

<Po.=^oi±^^^ln^^, (33) 

where the subscripts 1 and 2 refer to the corresponding examples. This is because the theory has a symmetry 
which is lost in the limit. 

Notice that the scalar field potential H29|) is well-defined at — 1 < fc < 1 and the model possesses the exact 
solution H3U|I which corresponds to the cosmological evolution (|28|l . The dynamics of the universe within the 
tachyonic models (|24l) and H31|l will be studied in detail in the next sections for general initial conditions. We 
shall show that both models can be extended to the interval < fc < 1, actually making them richer than the 
corresponding scalar minimally coupled theories. 
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3. In the third example we consider a perfect fluid whose equation of state is as follows: 

£+p = je^, 0<7<1, A > 1. (34) 

In this case the Hubble variable is 



( ^7 (2A-l)i) . (35) 



The scalar potential has the form 



,' 3V7(A~l)(^-V^o) \ ^ 7 / 3V7(A-l)(y-^o) ' ^^^^ 



The exact solutions are 



The tachyon potential is 



l-A 
1-2A 



^(i) = ±V7^^(^7(2A-l)i) +^0- (37) 



V{T) 

with the following exact solutions 



T{t) = ±^^—^-^(2X-l)tj +To. (39) 

The scalar field model with the potential (|36|l was studied in ■37^ while the corresponding tachyon model with 
the potential was discussed in (13^.14|. 



4. In our last example we consider the Chaplygin gas, described by the following equation of state: 

A 



P 



-, A>0. 



In this example the evolution h{t) is given only implicitly |43 | by the formula: 



1 



6^1/4 



In 



2 arctan 



^1/4 



(40) 



(41) 



However, as shown in (42||, the scalar potential can be reconstructed using the known explicit dependence of e 
on a and one gets 



= ^cosh3((^ - ifo) + 

The corresponding field configuration is also given implicitly: 

/ 



1 



cosh 3 ((/3 — ifo) 



(42) 



ip{t) = =F-arccosh 



A 



- 1 



v 



J 



Similarly, using the dependence of e on a one can reconstruct the tachyon potential: 

V{T) = VA^ const. 



(43) 



(44) 



It is easy to see [13 that the tachyon model with a constant potential is exactly equivalent to the Chaplygin 
gas model. Indeed, in the case of a constant tachyon potential the relation between the tachyon energy density 



-V^(T) 



'A. The Chaplygin 



Q and the pressure © is just that of the Chaplygin gas H4l)|) . where p e 

gas cosmological model was introduced in and further developed in l43l l44l l45l l46j| and many other papers. 
Comparison with observational data has also been extensively performed |47| . 
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IV. TRANSGRESSING THE BOUNDARIES 

We now take a step back and consider the problem of finding a tachyonic field theory admitting the same cosmic 
evolution as the one produced by perfect fluid with equation of state p = ks, where now k > 0. As we have stated, 
it is impossible to reproduce this dynamics using Sen's tachyonic action One way out is to introduce a new field 
theory based on a Born-Infeld type action with Lagrangian 



L = W{T)Vt^ - 1. (45) 
In this new field theory the energy and pressure are given by 

W{T) 



_ 1 ' 

and 



(46) 



P = W{T)VT^ - I. (47) 

The pressure (if well defined) is now positive. On the other side, the equation of motion for this field has exactly the 
same form © that was derived by Sen's action. 

Following the procedure described in Section 2, now applied to the Lagrangian H45fl . one gets the following potential 
corresponding to the equation of state p = ks (with fc > 0): 

Wm^^ ^ ^ (48) 

^ ' 9(i + fc)r2 ^ ' 

The exact solution of the field equations that reproduces the dynamics of the perfect fiuid is 



T{t) = VT+k t (49) 

(we restrict our attention to the region of the phase space where T > and T > 1). There are two other obvious 
solutions for this model, corresponding to other choices of the initial conditions; they also give rise to linearly growing 



fields: T{t) = t and T{t) = .y/l + i t. 

We would like to point out an interesting fact: the Lagrangian (|45|l which, together with the explicit form H48|) of 
the potential, describes the field theory corresponding to a positive value of k, is actually the same as Sen's Lagrangian 
(PI with potential (|24|) itself considered for positive k. Indeed, it is true on the one hand that in this case the potential 
(|24|l becomes imaginary. However, this can be compensated by considering the kinetic term in the region 1 — < 
so that the action as whole remains real. It can be re-interpreted as the product of two real terms: 

L = -V{T)^Jl -f^ = {^^fV{T)\/l - T2 = W{T)\/f^ - 1. (50) 

This model is introduced here as a pedagogical introduction to the model with the potential (|31|) . which we discuss 
in detail in the following section. The properties of the present model can be recovered in the limit A ^ and we 
will not further comment on it. 

V. DYNAMICS OF THE TOY TACHYONIC MODEL 

We now provide the analysis of the dynamics of the tachyonic model based on the potential (|31|l . In this case eq. 
© is equivalent to the following system of two first-order differential equations 

T = s, (51) 

s = -3VV{l-s'f/\s-{l-s^)yf, (52) 

where using the Friedmann equation Q we have expressed the Hubble variable h as a function of the variables T and 
s. The model has the following two exact solutions (we take Tq = without loss of generality): 

2 . , 3(1 + k)VAt , , 

ri(^) = — , arctansmh -^-^ , < i < oo, 53) 

3v/A(l + k) 2 

2 ( , 3(l + fc)VAt\ , , 

t = , TT - arctansinh -i — Q<t<co. 54) 

3VA(m)V 2 y 
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By inserting H31() into (|52() and by eliminating the time we obtain an equation for the phase space trajectories s = s{T): 
ds _ 3(1 ~s')^ (l-ik + 1) cos^ VA(i+fc)T \ 



dT . 3VA(l + fc)T I 1 - s2 

cm i 1 . \ 



sm 2 



3v/A(l + A;) 1 - ( 3 VA(1 + fc)r ^ [k + 1) cos^ ^^(1+^)^ + {k-l) 

cot I I , : ■ (55) 

/ 1 - (fc + 1) C0S2 2VAp)T 



2 s \ 2 



2 

In the phase plane (T, s) the solutions H53|) (|54|l correspond to arcs of the curve a (see Figs. ^ 121 and 0)): 



3v/A(l + /c)T 



s = vT+fc cos — ^ (56) 

The behavior of the cosmological radius for both solutions H53|l and (|54|l is the following: 

/ , 3VA(l + fc)t\™ 
a{t) = aa smh ^- ^ 1 . (57) 

As expected, we get back the cosmological evolution (|28|l determined by the equation of state (|27|l which was the 
starting point for constructing the potential (|31|l . To study the cosmological evolutions corresponding to all possible 
initial conditions we need to distinguish two different cases. When the parameter fc < there are no surprises and the 
associated cosmology is essentially driven by that of the exact solutions while when k is positive the model seems at 
first to be ill-defined. We will have again to go beyond the model itself and " transgress the boundaries" to see what 
its possible meaning can be. 

In the first case when — 1 < A; < the potential (|31|l is well-defined for 

277 

0<r< — ; (58) 

3v/A(l + fc) 

while the dynamics guarantees that 

-1< s < 1. (59) 

The system has only one critical point, namely 

so = 0, To = J" (60) 
3v/A(l + A:) 

The eigenvalues of the linearized system in the neighborhood of this point are 

Ai.2 = -^VA(l±fc). (61) 

Both of them are real and negative. Thus, this special point is an attractive node. It corresponds to a de Sitter 
expansion with a Hubble parameter 



Ha - VA. (62) 

The set of integral curves of 155|1 is symmetric under refiection with respect to the critical point H6U|I : any given 
integral curve and its node (n)-refiected one describe the same cosmological evolution. The curve a corresponds to 
the eigenvalue Ai = — |-\/A(l + fc), whose absolute value is the smallest of the two. It acts as a separatrix for the 
integral curves. Almost all curves approaching the node (|60|) end up there with the same tangent as a. The only 
exception is a second separatrix a' which corresponds to the eigenvalue A2 = — |\/A(1 — fc). The curve a' separates 
the bundle /' of the curves which do not intersect the axis s = from the bundle / of those which intersect it (see 
Fig. 1). The boundary of the rectangle defined by Eqs. H58|) and (|59|l describes a cosmological singularity. Indeed, 
the scalar curvature for a flat Friedmann universe is 



R = &{h + 2h^). 



(63) 
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Since from Eqs. (P), Ql and 101 one has that 

3 



h = -^h'^s^ (64) 



by substituting into (|63|l it foUows that 



i? = 3/.^(4-3.^) = ^^:™^. (65) 
Vl — 

Thus the scalar curvature R tends to infinity when approaching the boundary of the rectangle (with the exception of 
the corners, that will be treated separately). 

All integral curves end up in the node; let us see how they behave close to the boundary and begin with a right 
neighborhood of T = 0. There equation H55|l takes the form 

2{1 -s^)f^ (^fc)i/4, ^_Fis) ^gg^ 



dT sT \ 7T+-fc(l-s2)i/4y T 



If s < vT+fc, then s +oo as T — > 0. Therefore, the integral curves at s < vT+fc which get close to the T = 
axis rise almost vertically, climbing leftwards for s < and rightwards for s > until they get close to a whenceforth 
they reach a maximum and thereafter approach the de Sitter node (see Fig. 1). These are the curves of bundle 
/. If y/TTk < s <1, then s —oo as T ^ 0. Therefore, the corresponding curves at s > \/l + k which get close 
to the T — Q axis drop almost vertically until they get close to a whenceforth they also approach the de Sitter node. 
These are the curves of bundle /'. On the separatrix a one attains the point (0, Vl + fc) where s — 0. Symmetric 
considerations apply to the n-reflected curves, i.e. those which lie to the right of the separatrix a' . 
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But where do all the curves originate from? We first show that apart from (T, none of them can touch any point of 
the s-axis. Indeed, let us consider a point (T, s) close to the s-axis. Eq. (|66|) can be integrated backwards to give 

f dx 

T(s)=ro(so)exp / -— . (67) 



So 



Fix)- 



This equation shows that, if F{so) ^ 0, it is impossible to realize the condition To(so) = and therefore touch the 
s = axis on a given trajectory. The roots of the equation -F(so) = are s = yl + fc, ±1. A closely similar reasoning 
excludes the point (0, 1). The point (0, —1) can also be excluded since we should have T > in a neighborhood of 
such point; but s is negative and this contradicts the equation T = s. We are therefore left with the point (0, ^/l + k), 
where the exact solution a originates. 

Let us examine now the upper boundary of Fig. ^ In a small neighborhood of the point (T,, 1) (where is in the 
domain (jSHl) Eq. ifK^ can be replaced by the following approximate equation: 



ds 



3(2(1 - ,s))3/4VA (^1 - (1 + k) cos2 VMi+fe)T. ^ 



dT . 3jA{l + k)T, 



(68) 



sm- 2 



This equation is not Lipschtzian. The upper integral is the trivial solution s = 1 while the lower integral is 

= 1 for T < 

i 1 - C{T,){T - T,)4 for T > T, ^""^^ 
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2 M n , ^^ „„„2 3\/A(i+fc)r. 



where 

'"^■^*^ ~ 32 ~^JJW+Wr. ' ^"^^^ 

sm 

The intermediate solutions stay constant at s = 1 for a while and then leave the s = 1 line at a value T** > T^. 
Therefore, from each point (T»,l) there originates only one integral curve behaving as H69|) . In particular, the 
separatrix a' originates at a point {T„i , 1) (the value of T^i being unknown). 
The condition of cosmic acceleration is 

— = -(e + 3p)>0. (71) 
a 

For the tachyon cosmological model, using formulas (@J) and lO this condition can be reexpressed as 

.2 < I (72) 

Therefore, if fc < — i all cosmological evolutions undergo an initial phase of deceleration followed by an accelerating 
one. On the other hand, if fc > — i all evolutions whose originating point T, is larger than a critical value (which 
depends on k) have two epochs of deceleration and two epochs of acceleration. This happens in particular in the 
limiting case when fc = 0, as shown in Fig. |21 

We consider now the case < fc < 1. Eq. (|31|l shows that the potential V{T) is well defined only in the interval 
(T3,r4) where 

2 1 2 / 1 \ 

T3 = , arccos T4 = , tt — arccos , . (73) 

3^/(1 + fc)A 3^(l + fc)A V VlTk) 

Together with the condition — 1 < s < 1, Eq. 1)73(1 defines the rectangle where we study the model at first. There 
are now three fixed points. One of them is the attractive de Sitter node H60() whereas the other two are saddles 
corresponding to the maxima of the potential at coordinates 



2 /I — fc 2 / 1 — k\ 

Ti — , arccos \ / , T2 — , n — arccos \ / . (74) 

3y/{l + k)A Vl + fc 3^/(TTk)A\ \l + kj ^ ^ 

(s=0). They give rise to an unstable de Sitter regime with Hubble parameter 



We first analyze the behavior of the trajectories in the vicinity of the line T — T3 and set T = T3 + T, with T small 
and positive and s 7^ ±1. With these conditions the model is described by the approximate equation 

s« ^, 76 

2T 

which implies that the trajectories passing close to the boundary T — drop steeply down without crossing it. The 
"physical" reason for this behavior is the vanishing of the potential a.t T — (to our knowledge this is a novel feature 
of our model) ; indeed, the structure of the tachyonic action implies that the " force" is proportional to the logarithmic 
derivative of the potential and this is infinite a,t T = T3. The impossibility of crossing T ^ remains true also 
without coupling the tachyon to gravity. On the other hand, in contrast with the situation encountered before, the 
geometry is regular at T = T3: the vertical boundaries of the rectangle are not curvature singularities, because the 
potential does not diverge there. Actually, the curvature scalar R vanishes there because of the vanishing of the 
potential. Instead, the horizontal sides are still singular. The situation is exactly as before and there is one integral 
curve which originates at {T^,, 1) whose behavior is again given by (|(j9|l . H7U|) . Now, however, due to the positivity of 
the parameter fc, the nonnegative function Cfe(T*) vanishes at T = Ta and at T = T4, is maximal at T = Ti and at 
T = T2, and has a positive local minimum at T = Tq. 
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FIG. 3: Potential for fc > (here k = 0.44). 



Now consider the behavior of the trajectories in the upper left vertex of the rectangle P ~ (Ts, 1). Setting s = 1 — s 
and T = + T, with T and s small but not zero, from Eq. H55|l we get the approximate equation 



dT T 

where A = 3^/^23/4 A5/8(l ^ kf/»k-^/^_ The general solution is 

—{AT + Bff, 
256^ ' ' 

where B is an arbitrary constant. When B the leading behavior is 

S = DT, D>0. 



(77) 



(78) 



(79) 



This means that there are trajectories which gush out from the point P in all possible directions, except of course 
vertically (see Fig. 4). If B = Eq. (JTH)) becomes 



~s = — f 5. 
256 



(80) 



This equation describes the leading behavior of a curve p which separates the trajectories of type H79() from those 
originating at the points (T*, 1). We point out that P is a cosmological singularity for the curve p (which therefore 
originates at P), whereas it is regular for the curves H79|) . Indeed, we see from Eq. ^ that when T ^ behaves 
as l/T^ along p, while it is finite along trajectories H79|) . 

Now consider the upper right vertex Q — (T4, 1). By setting as before s = 1 — ,5 and T = T4 — T we get the 
approximate equation 



dT T 



where A is the same as above. The general solution is 

1 



s ^ —{-AT + BfT, 
256^ ^ > ' 



(81) 



(82) 



where B is an arbitrary positive constant. Therefore the trajectories enter point Q along all possible directions except 
vertically and horizontally. We now classify the behavior of the trajectories in the interior of the rectangle. First note 
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FIG. 4: Phase portrait evolution for A: > (fe = 0.44). 

that there are five distinguished trajectories (separatrices) : a which connects P with the node; t which connects P 
with the saddle (Ti, 0); ^, which is the curve which enters the saddle (r2, 0) from above (We are unable to say whether 
^ originates at some point (T*, 1) or it belongs to the family H79|) . or it coincides with the curve p defined by equation 
(|80|l : for this reason we have not tried to draw curve p in Fig 01; t/;, which originates from the saddle (Ti,0) and 
enters the node with tangent defined by the eigenvalue A2 — — |\/A(1 — k) (see (j^ l: finally, x which connects (T2, 0) 
with Q (see Fig. 0J. Each of these separatrices has its own n-rcflected counterpart (denoted by the same symbol). 

Corresponding to the separatrices cr, r, ^, ip and x can distinguish four bundles of qualitatively different trajec- 
tories. 1) The bundle of trajectories limited by r, ip and a: they originate from P and enter the node along ij) and 
from above: 2) the bundle /^^ of trajectories limited by a, ^ and ip : they enter the node along ip from below; 3) the 
bundle f^^^ of trajectories limited by ^, x and the horizontal line s = 1; they stream into Q\ 4) the bundle /^^ of the 
curves which are limited by r, x a-nd the vertical line T = T^: they gush out of P and stream into Q' . 

Now we have to face a problem that we have not mentioned so far. The question is the following: it takes a finite 
proper time for the fields (and the universe) to get from anywhere to the corner Q or Q' on a trajectory of bundle /^^^ 
or of bundle /^^. But these corners are not critical points of the dynamical system and, furthermore, the universe 
does not experience any singularity by getting there along these trajectories. Similar remarks apply to the past of 
the trajectories originating from P and P' . If the model could not be extended to follow these trajectories outside 
the rectangle where it has been originally defined it would be useless. However, we now show that this extension 
is actually possible. Indeed, one can see by inspection that the field equations are well defined also inside the four 
semi- infinite strips defined by the following inequalities: < T < Ts, with s > 1 01 s < —1; < T < — . •^^ with 

t^{l+k) 

s>lors<— 1. Then, following the strategy sketched in Section ITVI we introduce a "new" Lagrangian: 



L = W{T)Vt^ - 1, 



(83) 
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where the "new" potential is given by 




(84) 



The "new" Lagrangian comes itself out of Sen's action (Eqs. (|^. 13111 ) by the previous trick: 



V{T) — > W{T) = iV{T) 



(85) 
(86) 



Vi -T2 — > iVi-t^ = Vt^ - 1 



i.e. both the "old" kinetic term and the potential become imaginary but their product remains real. It follows 
from the dynamics (|55|l that the expressions under the square root in the potential and the kinetic term change sign 
simultaneously. In other words, in the phase diagram no trajectory can cross any side of the rectangle. Thus, all the 
quantities characterizing the model keep real during the dynamical evolution. White regions in the phase diagram, 
where the Lagrangian and other quantities would become imaginary, are forbidden! The product H83|l , which amounts 
exactly to the "old" Lagrangian, can be interpreted in terms of the "new" kinetic term and potential that are both 
real. This makes it clear why the "new" Lagrangian gives rise to field equations in the above four strips, which are 
the same as in the interior of the rectangle. At first glance one might have the impression that there is a freedom of 
choice of sign for the new Lagrangian H83(l . or, in other words, that one may choose opposite signs in Eqs. (|86|l . JHSJ- 
However, this is not so: the choice of sign in Eq. H83|l is determined by the requirement of continuity of the Einstein 
(Friedmann) equations when passing from the rectangle to the stripes. As anticipated in Sec. IIVI both the energy 
density and the pressure are positive in the considered strips. 

However, there is an important difference between the present situation and the one described in Sec. IIVI Here we 
are dealing with one single model (fc is fixed). One is compelled to make the extension of the model and, contrary to 
the A = case, the two different "phases" (i.e. negative and positive pressure) are found within the same model, at 
different stages of the cosmic evolution, one phase in the rectangle, the other in the strips. In the following we give 
the precise mathematical meaning of this extension. This opens the way to the study of a new class of tachyon field 
theories. 

We start by describing the behavior of the trajectories in the lower left strip (see Fig. 0J. Since T < — 1, the 
evolution along any given trajectory will lead us in a finite amount of time to either hit the vertical line T = or 
to approach a vertical asymptote T — Tg, with < Tg < T^. The first alternative does not take place, whereas all 
values of Tb in the indicated range are allowed, with the exception of Tb = 0. Indeed, in the vicinity of T = Eq. 
(|55|l takes the following form 



Now, assume there is an sq < — 1 such that lim^^o •s(r) = sq. Then, the analogue of Eq. (|67|) gives a contradiction: 
the function F{s) cannot vanish since the r.h.s. of Eq. 187|l is positive. 

The leading behavior of the solutions of Eq. (|S7|l for T ^ 0, s ^ — oo, is given by 



where Tb is an arbitrary positive constant. It follows that the line T = is not an asymptote for any trajectory but 
of course there are trajectories whose asymptote is as close as we like to the line T = 0. 

Therefore, the trajectories inside the considered strip can be parameterized by the value Tb of the coordinate of 
their vertical asymptote. We now discuss how they behave in the neighborhood of s = — 1. As before we can divide 
the trajectories inside the strip into three families according to their leading behavior 



ds _ 2(s2 - l)3/4fcl/4 2(.s2 - 1) 



(87) 



dT VTTkT sT 




1 



(88) 



s 



1 + E{T~T3), E>0, 



(89) 



s 



1 + (T-Ta)^ 

256^ ' ' 



(90) 
(91) 



s 



1 + C(T^){T - nf, < < Ta 



Trajectories of type H89|l fan out from point Q' into the strip at all possible angles. Trajectories of type (|91f) can be 
thought of as originating from (T*,— 1). Note that, just as before, the function C(T,) approaches (minus) infinity as 



15 



— > 0. It approaches zero as T» —^ T3. Curve H90(l separates the famihes (|89|l and 1)91(1 . The coordinate Tb of the 
asymptote is an increasing function of . As — + it attains a value characterizing the asymptote of curve H9U|I , 
beyond which it becomes an increasing function of E. 

Let us go back to physics and consider the behavior of the cosmological radius a{t) when the tachyon field T{t) 
tends to Tb along the solution of the field equations. The Friedmann equation implies that /i^ — > and that h — > —00 
(see Eq. (|64ll ) . Therefore, the scalar curvature H63|l diverges and the universe reaches a cosmological singularity in a 
finite time. 

This is an unusual type of singularity which we call big brake. Indeed, since a/a = /i + /i^, in a big brake we have 
that 

a —00, 

a — > as < 00 (92) 

In other words, the evolution of the universe comes to a screeching halt in a finite amount of time and its ultimate 
scale depends on the final value Tb of the tachyon field. 

We now turn to the behavior of the trajectories in the upper left strip. In the vicinity of T = the equation for 
the trajectories takes the form (|87ll . The coefficient of T^^ vanishes at the values 



1 



5=1, VT+k, + (93) 

As explained earlier it is only at these values of s that the trajectories can leave the line T = 0. In addition, an 
analysis similar to the one performed earlier shows that the only trajectory starting at (0, 1) is the line s = 1, and 

that the only trajectory starting at (0, y/1 + k) is the curve a. All other trajectories start from ^0, + ■ Eq. 
(|87|l shows that as T the derivative ^ approaches —00 for \/l + k < s < W whereas it approaches +00 



for 1 < s < \/l + k and s > 1/ (see Fig. 0J|. To study the behavior in the neighborhood of 0, \/ ) we set 



s = y ^-Y^ + s with s small. In the approximate equation for s{T) it is necessary, besides the leading term , to keep 
the term proportional to T. Thus we get 



ds 2(1 -fc) S 9A 1 + k fk + 3 



dT 1 + k T 8 
The general solution of this equation is 

\3/2 



T. (94) 



9A(l + fc)3/2(/c + 3)^4. \ 



S= D- ' LT— T^, (95) 



32A:5/2 

where D is an arbitrary constant. If Z? 7^ we have the leading behavior 



1 2(1-*=) 

1 + - + DT^TT (96) 
k 



Therefore, if fc < ^ the trajectories start from ^0, y 1 + with horizontal tangent, whereas they start with vertical 
tangent when fc > i. If fc = they are born with any possible tangent. When Z) = we must go to the next order: 



1 _ 9A(l + fc)3/^(fc + 3) , 

32fc'V2 ^ ■ ^""i 

This curve acts as a separatrix for the curves having positive and negative D respectively. Regarding the curves of 
the strip which lie below cr they depart from (T*, 1) and behave in the neighborhood of this point as 

s = l-C7(^,)(T-T0^ 0<T, <r3, (98) 

where the function C(T*) is the same as in H91f) . Regarding the behavior of the trajectories in the neighborhood of 
T = T3 a simple analysis shows that they stream into P at all possible angles (except vertically and horizontally). 
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These properties show that each curve of type (|96|l with positive D attains a maximum somewhere between T = 
and T = whose height is an increasing function of D which tends to infinity as D +00. 

So far, we have analyzed the behavior of the trajectories in two distinct regions: the rectangle and the four strips. 
Now, it has to be noted that the trajectories in the rectangle which leave P (or P') at all possible angles in the open 
interval (0, 7r/2) and the trajectories which enter Q' (or Q), again at all possible angles, are incomplete, since the 
vertices of the rectangle are not cosmological singularities for these curves. The same is true for the trajectories in 
the strips which enter P (or P') and leave Q (or Q'). This circumstance, and the fact that the equation of motion 
for T is the same in the rectangle and in the strips, indicates that it must be possible to extend and complete the 
above set of trajectories by continuation through the vertices of the rectangle. Precisely, the trajectories entering P 
(or P') from the upper left (or lower right) strip shall be continued into the trajectories entering into the rectangle 
from P (or P'). Similar remarks apply to the corners Q and Q' . The uniqueness of this continuation procedure can 
be proved by applying the cr-process method of resolution of singularities (see, e.g. ^3). This amounts to blowing 
up (unfolding) the vertices of the rectangle by transition to suitable projective coordinates, the use of which removes 
the degeneracy of the vector field at these points. We do not give the mathematical details. 



VI. COSMOLOGY 



The rich mathematical structure that we have exhibited in the previous section gives rise to the possibility that 
cosmologies that have very different features coexist within the same model. In this sense tachyonic models are richer 
than the "corresponding" standard scalar field models. 

We now characterize all the evolutions in our model that are portrayed in Fig. 0] We start from the trajectories 

originating from the cosmological singularity ^0, y^l + and which are characterized in the neighborhood of the 

latter by formula (|95|l with the parameter D positive and very large. One such trajectory, as soon as it leaves the 
singularity, raises steeply upwards until it attains some maximum value. Henceforth it turns steeply downwards, 
enters the rectangle at P at some small angle a with the vertical axis T = T3 and moves towards Q' . Upon reaching 
this point it enters the lower left strip and eventually ends up in a finite time tg in a big brake corresponding to some 
field value Tb very close to T3 (inf ts =0). 

As D decreases the height of the maximum decreases accordingly, a and ts increase and Tb decreases. Eventually, 
at some critical value Dc of the parameter D, the trajectory degenerates inside the rectangle with the separatrices r 
and X entering and, respectively, leaving the saddle point (Ti, 0). The curves for which D > Dc belong to the bundle 
/^^. For D < Dc we get the bundle f^. These trajectories correspond to evolutions which are asymptotically de 
Sitter with a Hubble parameter Hq = VA- The upper bound of D for the curves of bundle is Dc and corresponds 
to the separatrix t and ip. This means that a tiny difference in the initial conditions will end up into dramatically 
different evolutions: one universe goes into an accelerating expansion of the de Sitter type and the other ends with a 
big brake. This should not be confused with a chaotic behavior: the two evolutions are almost indistinguishable for 
a very long time and then suddenly diverge from each other. As D approaches —00 we end up into the separatrix a 
originating at the cosmological singularity (0, Vl + k). 

Proceeding further we encounter those evolutions which start at (T*,l), < T^. In the neighborhood of the 
starting point they behave according to Eq. (j^Hl), then enter the rectangle through P passing above the curve a 
from below. In the limit = T3 we obtain the trajectory p. As increases further beyond T3 the curves detach 
themselves from the axis s = 1 according to the equation H98|l now with < < T4. At some critical value 
of T*,0 < < T4, the trajectory degenerates into the separatrices ^,1/1 and x- The curves for which < T* < 
form the family /^^ and they asymptotically arrive at the de Sitter node (To, 0) from below the axis s — 0. Those 
for which < T^, < T4 form the family /^^^; they enter the upper right strip through Q and eventually end up in a 
finite time in a big brake. Finally, the trajectories that detach themselves from the axis s = 1 at T, > T4 bend again 
upwards and they form another bundle that we denote by They also end up in a big brake. 

The evolutions corresponding to the trajectories of the different bundles can also be studied in terms of the qual- 
itative behavior of the Hubble parameter h{t) as a function of time (see Fig. [SJ. In our model, using equations 
(|ST|l and lIHSI) we get 



1 



1-S2(t) 



(99) 



Since h = —^h^s^, h{t) is positive and strictly decreasing, except at those values of t for which s{t) = 0. 




Close to the initial singularity, the trajectories can be divided into two classes depe nding on th e singular point from 
which each of them starts: the class (A) of trajectories originating at the point (0, ^^(1 + k)/k) and the class (B) of 
trajectories which start from the points (T,, 1). The two classes are separated by the curve a. The leading behavior 
of hit) in the neighborhood of the singularity at t = for the trajectories of class (A) is: 

= 37T?feTi-^(¥)'^*^+---' for i?^0; (100) 



3(l + fe) I 

3(l + fe) I 32fc 

For trajectories of class (B) the leading behavior is universal (it does not depend on k) 



hit)^jrtm\-^-^M+..., for D^O. (101) 



/i(t) = | + ... (102) 

To study the behavior of h{t) at later times, one needs to examine each bundle / separately, which, in turn, requires 
first plotting h{t) for the different separatrices cr, r, and x- 
We have from Eq. (|57|l that h{t) for the separatrix a is given by 



Kit) = VAcoth (^^^ 



k)t 



(103) 



In particular, in the neighborhood of t = 0, 



, 2 1 A(l + /fc) , , 

Mt)^^(^^ + ^^ + ... (104) 

to be compared with Eqs. pOO|) . (|101|l and (|102|) . Furthermore, ha-{t) is strictly decreasing and, as i — *■ oo, it 
approaches a stable de Sitter expansion with Hubble parameter \/A- 
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Regarding the unstable cosmological evolution hr {t) , for small t we have 

, 2fc 1 Dck f l + k\^^ i-3k 

and limf^oo hr{t) — Hi (see Fig. Eland Eq. H75(l ^. Qualitatively, hr{t) behaves similarly to hcr{t) but the asymptotic 
value of the Hubble parameter is higher. As for the separatrix ip, the corresponding cosmological evolution is unstable 
and non singular: the Hubble parameter decreases steadily from an unstable de Sitter regime h^{—oo) — Hi at 
large negative times to a stable one h^(po) = \/A at large positive times. /i^(t) and /ir(i) have qualitatively similar 
behaviors at large times while they are different at small times. Finally, h^{t) decreases steadily from an unstable 
asymptotic de Sitter value h^{—oo) — Hi at large negative times to the final big brake singularity h^(tB) = 0, with 
h^{tB) = — oo, for some suitable ts- 

For a trajectory of bundle which lies close to the sepratrices t and tp, h{t) behaves according to H1UU|) close to 
the initial singularity, with D < Dc, {Dc — D) small. Then h{t) decreases steadily and remains close to hr{t) for a 
very long time during which s (and hence h) gets close to zero. Therefore, in this regime the evolution simulates an 
asymptotic de Sitter expansion with Hubble parameter close to the value Hi. Eventually, however, s starts increasing 
again and the graph of h{t) bends downwards away from hr{t) and approaches asymptotically the stable de Sitter 
regime with Hubble parameter VA. Instead, for a trajectory of bundle which lies very close a the behavior Ij 100(1 
of h[t) at small times is characterized by a value of the constant D which is negative and very large; the graph of 
h{t) parts only slightly from the graph of hcr{t) and the asymptotic value a/A is approached much earlier. Other 
elements of display behaviors which are intermediate between those described above. The time dependence of h 
for the trajectories of bundle /^^ is qualitatively similar to the one relative to the curves of the differences being 
the following: the small time behavior is given by (|102|) . and for each trajectory of /^^ there is a value Iq of t (which 
depends on the particular trajectory) for which s(io) = so that /i(to) — 0- 

Now consider a curve of bundle /^^^ which lies very close to the separatrices f and x- For such a curve h{t) remains 
very close to h^{t) for a very long time, s{t) decreases getting close to zero and the evolution simulates again an 
asymptotic de Sitter expansion with Hubble parameter. However, eventually s{t) starts increasing indefinitely and in 
a finite (though long) time the cosmology end up in a big brake. Moving to curves that are farther and farther away 
from ^ and x, the value of T at the initial singularity moves to the right towards the value 27r/(3-\/A(l + k)), h(t) 
decreases more and more steeply and the big brake time is tends to zero (when gets larger than T4 the trajectories 
switch from bundle /^^^ to bundle /^). 

Like those of bundles /^^^ and the evolutions of bundle f^^ give likewise rise to a big brake and the behavior 
of h{t) is qualitatively similar in the two cases. However, contrary to what happens for bundles /^^^ and for each 
trajectory of bundle /^^ there is a time t = (which depends on the particular trajectory and spans the whole open 
half line) for which s{to) = and hence /i(to) = 0. Therefore, even though the big brake time can get as close to zero 
as one likes for the curves of type IV (for such curves, is is a decreasing function of the parameter D in Eq. ((l()()|l . 
with Imiu^u^ tB^D) = 00 and limu^oo tsiD) — 0) and hence h{t) may decrease to zero very steeply, there is always 
some intermediate time at which h{t) is momentarily stationary. Remark that, as one can see from Fig. 5, there is 
nothing peculiar in the behavior of h(t) at the times when vertices of the rectangle are crossed. 

The possibility of cosmological singularities characterized by the divergence of the second time derivative of the 
cosmological scale factor has also been considered in in the brane cosmology context. 

We conclude this section by another simple example of cosmology sharing this property. Let us consider the flat 
Friedmann universe filled with a perfect fluid with a state equation similar to ()40|l : 

A 

p=- (106) 

where A is positive. We can call this fluid the "anti-Chaplygin gas". This equation of state arises in the study of the 
so-called wiggly strings [sol Isif . The dependence of the energy density on the cosmological radius is given by 



^-A, (107) 



where B is a positive constant. At the beginning of the cosmological evolution e '-^ yBjar , like in the dust-dominated 
case. Now there is a maximal value possible for the cosmological scale 

5X1/6 



aj. = ( - ) . (108) 
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that is attained in a finite cosmic time tp. The behavior of a(t) in the vicinity of the maximum is the following 

a{t) ^ap- C{tF - C = 2-^/335/3(^5)1/6^ (^Og) 

Since d{tp) — while a{tp) — —00 we are back into a big brake cosmological singularity. 

Thus, a big brake singularity can be found in an elementary cosmological model (though based on an exotic fluid). 
The difference with the tachyonic model is that in the latter there are evolutions culminating in a big brake which co- 
exist with other evolutions giving rise to an infinite accelerated expansion. These two types of evolutions correspond 
to different classes of initial conditions. 

It maybe worth mentioning that also for A; > one can construct a scalar field model; indeed the potential displayed 
in Eq. (|29|l is not restricted to fc < 0. But this model has a much poorer dynamics: all the trajectories tend to the 
de Sitter attractive node point. 



VII. FINAL CONSIDERATIONS: FATE OF THE UNIVERSE 

The study of the fate or the future of the universe is rather popular nowadays (53~Iz3j B,nd, as was predicted more 
than 20 years ago by F.J. Dyson "s^l, "eschatology" has now become part of the cosmological studies. This sort of 
studies include of course also biological aspects of and the fate of consciousness in the different cosmological scenarios 
(see, for example, IHH I8(tI | 1 but our goal is much more modest and we shall concentrate on the purely geometrical 
facet of the topic. 

There are three mainly studied possible scenarios for the future of the universe, intensively discussed in literature: 
an infinite expansion; an expansion followed by a contraction ending in a "big crunch" ; an infinitely bouncing or 
recycling universe. 

The present set of observational data seems to favor the first scenario: the data are quite compatible with a flat 
universe with a positive cosmological constant. 

On the other hand a negative value of the cosmological constant fits better with string theory (see, e.g. [l^, ) ■ 
The presence of a small negative cosmological constant could be made compatible with the present day cosmic 
acceleration, provided there are some fields or types of matter responsible for this acceleration. However, in this 
context the cosmological radius sooner or later will start decreasing and the expansion will be followed by a contraction 
which would normally lead the universe to a cosmological singularity of the big crunch type. 

The third scenario of an infinitely bouncing or recycling universe appears to many to be very attractive because it 
opens an opportunity to escape the "frying" in the big crunch and the "freezing" in the infinite expansion case. In 
this scenario, during the process of contraction there will be two opportunities: collapse or bounce with subsequent 
expansion. The choice of one of these opportunities depends on the initial conditions and this dependence has usually 
a chaotic character [Til Fl Fljj^ iTTj . 

Furthermore, one can show 63] that for every model of dark energy describing an eternally expanding universe one 
can construct many closely related models which describe the present stage of acceleration of the universe followed by 
its global collapse. However, these models corresponding to eternally expanding and collapsing universes are different 
and have different values of their basic parameters. 

One interesting feature of our toy tachyon model is that the first and the second scenario coexist in its context. 
Depending on initial conditions, some correspond to eternal expansions of the universe which approach asymptotically 
a pure de Sitter regime, while others end their evolution at the cosmological singularity. 

The second distinguishing feature of this model is that the singularity at which the universe ends its evolution is 
not a standard big crunch singularity. Instead, it corresponds to a finite non zero cosmological radius at which the 
Hubble parameter is finite and the deceleration parameter is infinite and has positive sign. We have called this fate 
"big brake". The prospect of hitting the cosmological singularity during expansion has been also discussed in Ref. 
|57j . where the singularity corresponds to an infinite value of the Hubble variable and the cosmological radius. This 
scenario is known as "big rip" or "phantom cosmology" (see, e.g. {65, 67, 68, 69] ). 

The third distinguishing aspect of our model is the fact that the regions of the phase space corresponding to different 
types of trajectories are well separated and the dependence of the cosmology on the choice of initial conditions is quite 
regular (non chaotic). A remark is in order here: the chaoticity of the classical dynamics hinders the application of the 
WKB approximation and, hence, undermines the basis of the majority of results of quantum cosmology In our 
model quantum- cosmological schemes of the traditional type |8ll js^, |83, ImI can be attempted. The corresponding 
wave function should describe a probability distribution over different initial conditions for the classical evolution of 
the universe. The quantum evolution of the universe in our toy model might be expressed in the language of the many- 
worlds interpretation of quantum mechanics [85l l86j | . In this framework one can say that the wave function of the 
universe describes the quantum birth of the universe and subsequently the process of the so called "classicalization" 
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(see e.g. Ref. |83). During this process the wave function sphts into different branches or Everett worlds corresponding 
to different possible classical histories of the universe. The peculiarity of our toy model consists in the fact that some 
of these branches describe eternally expanding universes while other branches correspond to universes which end their 
evolution hitting the cosmological singularity. 
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